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SUMMARY 


A  semiempirical  profile-correction  factor  is  discussed  which  enables 
the  estimation  of  the  wave  drag  due  to  thickness  at  supersonic  speeds 
for  three-dimensional  wings  with  arbitrary  airfoil  section  (subject  to 
the  restrictions  of  thin-airfoil  theory)  through  use  of  previously 
calculated  drag  coefficients.  Application  of  the  proposed  correction 
factor  to  the  known  drag  coefficients  of  some  rhombic-profile  wings 
yielded  estimates  for  the  drag  coefficients  of  parabolic-arc-profile 
wings  that  were  in  good  agreement  with  theoretically  calculated  values. 

It  is  expected  that  satisfactory  estimates  can  be  obtained  for  many 
combinations  of  plan  form  and  profile  by  judicious  use  of  the  proposed 
profile  correction,  especially  at  speeds  for  which  the  wing  leading  edge 
is  supersonic. 


INTRODUCTION 


A  number  of  papers  dealing  with  the  linearized-theory  calculations 
of  supersonic  wave  drag  at  zero  lift  (that  is,  drag  due  to  thickness)  for 
three-dimensional  wings  have  been  published.  (See  references  1  to  12, 
for  example.)  As  a  result,  extensive  theoretical  data  are  available  for 
the  wave-drag  coefficient  at  supersonic  flight  speeds  for  wirig  plan  forms 
that  have  arbitrary  sweepback  (or  sweepforward) ,  taper  ratio,  and  aspect 
ratio.  The  profiles  of  these  wings,  however,  are  in  general  restricted 
to  symmetrical  double  wedges,  rhombuses  (symmetrical  double  wedges  with 
maximum  thickness  located  at  50  percent  chord),  and  biconvex  sections 
composed  of  two  symmetrical  parabolic  arcs.  Other  types  of  thin  small- 
slope  airfoil  sections  require  a  great  deal  of  mathematical  and  compu¬ 
tational  labor  and  have  been  avoided  in  the  theoretical  analyses. 
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Inasmuch  as  these  other  profiles  may  prove  to  be  of  interest 
because  of  structural  advantages ,  low-drag  characteristics,  and  other 
reasons,  it  would  be  desirable  to  be  able  to  obtain  estimates  of  the 
drag  for  these  cases  without  having  to  resort  to  further  theoretical 
investigations.  The  purpose  of  the  present  paper,  therefore,  is  to 
discuss  the  application  and  limitations  of  some  semiempirical  profile 
corrections  which  will  enable  the  estimation  of  the  wave-drag  coefficient 
for  three-dimensional  wings  with  various  profiles  by  utilizing  previously 
calculated  drag  coefficients. 


SYMBOLS 


M 

0 

A 

X 

A 

t/c 


free- stream  Mach  number 
Mach  number  parameter 


( \fw *  -  1 ) 

aspect  ratio  ^Wing  span)^/Wing  areay 
taper  ratio  (ratio  of  tip  chord  to  root  chord) 
sweep  angle  of  leading  edge,  degrees 

thickness  ratio  of  section  in  free-stream  direction  (maximum 
thickness  of  given  section  profile  divided  by  its  chord) 

wave-drag  coefficient  of  airfoil  section 
Drag  per  unit  span 


0 


Dynamic  pressure  X  Chord 


0 


wave-drag  coefficient  of  three-dimensional  wing 
Drag _ 


0 


5) 


Dynamic  pressure  X  Wing  area 
airfoil  section 

Use  of  subscripts  is  indicated  or  explained  in  text 


ANALYSIS  AND  DISCUSSION 

As  mentioned  in  the  introduction,  the  purpose  of  the  present  paper 
is  to  discuss  some  semiempirical  profile  corrections  which  will  enable 
the  estimation  of  the  wave-drag  coefficient  for  three-dimensional  wings 
with  small-3lope  profiles  by  utilizing  previously  calculated  drag  coeffi 
cients.  It  appears  that  profile-correction  factors  that  are  especially 
suitable  for  supersonic-leading-edge  conditions  can  be  obtained  from 
available  data  and  formulas  without  much  difficulty. 
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For  example,  suppose  the  drag  coefficient  is  desired  at  a 

free- stream  Mach  number  Mx  for  a  wing  with  aspect  ratio  Ax,  taper 
ratio  X1;  sweepback  or  sweepforward  of  leading  edge  A^,  section 
profile  sx,  and  a  constant  thickness  ratio  (t/c^  along  the  span.  If 
the  drag  coefficient  CDq  of  a  geometrically  similar  wing  at  the  Mach 
number  with  profile  sQ  is  known  (that  is,  if  the  drag  coefficient 

is  known  for  a  wing  with  all  characteristics  other  than  the  profile  s 
\e\_ca  those  of  bne  desired  wing),  then  the  following  expression 

LCDA)  x  Profile-correction  factor 

A1  ■ 

X1 
A1 

(V  c)x 

so 


where  the  subscripts  outside  the  parentheses  indicate  wing  characteristic 
and  Mach  number.  In  addition  to  the  leading-edge- sweepback  requirement, 
the  sweepback  of  the  maximum- thickness  line  for  the  two  wings  should  be 
as  nearly  the  same  as  available  calculations  permit.  Two  simple  correc¬ 
tion  factors  suggest  themselves;  these  are  stated  and  discussed  in  the 
following  paragraphs.  ' 


Consider  first  a  three-dimensional"  correction  factor  obtained 
by  forming  the  ratio  of  the  wave-drag  coefficients  of  two  wings  of 
equal  taper  ratio  (for  example,  untapered  wings  for  which  X  =  l)  with 
the  same  aspect  ratio,  sweepback,  and  thickness  ratio  and  at  the  same 

profilf1  r  &a  dhfhdeS+aed  Wing'  °ne  untaPered  wi-ng  having  the  desired 
p  01 lie  s1  and  the  other  untapered  wing  having  the  profile  sQ; 

that  is,  form  the  ratio  j  .  The  reason 

V  D0/M1,A1,Aq,  (t/c)1,X=l 

the  Xml  case  was  chosen  is  that  fairly  complete  theoretical  data 
exist  for  both  the  symmetrical  parabolic-arc-profile  and  the  symmetrical 
ou  le-wedge-profiie  untapered  wings.  This  approximation  was  suggested 
to  the  author  of  reference  9  and  was  applied  in  that  paper  to  compare 
tne  calculated  value  (based  on  a  theoretical  analysis)  for  the  parabolic 
arc -profile  tapered  wing  with  the  approximate  value  of  obtained 

by  using  previously  known  data  for  CD  (from  reference  6) 
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^rom  re£erence  8),  and  j (from  reference  5)  • 

Excellent  agreement  was  obtained  for  the  cases  considered.  This  factor 
can  be  considered  as  essentially  correcting  for  taper  ratio,  since  the 
untapered-wing  value  for  the  profile  under  consideration  must  be  known. 
By  rewriting  the  preceding  equation  with  the  correction  factor  inserted, 
this  point  can  be  clearly  illustrated: 


Inasmuch  as  only  the  symmetrical  double-wedge,  rhombic,  and  symmetrical 
parabolic-arc  sections  have  been  treated  in  detail  for  untapered  three- 
dimensional  wings,  and  since  even  less  detailed  calculations  are  avail¬ 
able  for  tapered  plan  forms,  the  use  of  this  type  of  correction  factor 
is  very  limited  and  restricted.  Also,  because  most  of  the  theoretical 
solutions  are  based  on  linearized  theory,  such  solutions  for  the  double¬ 
wedge-  and  rhombic-profile  wings  exhibit  an  additional  "drag  peak" 

(in  violation  of  the  small-perturbation  assumptions  of  the  linear  theory) 
whenever  the  Mach  number  is  such  that  the  Mach  lines  parallel  the 
maximum- thickness  ridge  line.  (This  condition  is  termed  a  sonic  line  of 
maximum  thickness  since  the  free-stream  flow  component  normal  to  that 
line  is  sonic.)  Hence,  good  agreement  with  experiment  or  with  more  exact 
calculations  would  not  be  expected  at  or  near  Mach  numbers  for  which 
this  condition  is  present.  It  might  be  added  that  equally  erroneous 
drag  peaks  result  from  using  linear  theory  whenever  a  line  connecting 
the  points  of  appreciable  wing-slope  discontinuities  becomes  sonic. 

For  example,  the  biconvex-,  rhombic-,  and  double-wedge-profile  wings  all 
exhibit  drag  peaks  at  Mach  numbers  for  which  the  leading  and  trailing 
edges  become  sonic;  hence  linearized-theory  calculations,  or  approxi¬ 
mations  involving  the  use  of  such  calculations,  at  or  near  these  Mach 
numbers  should  not  be  expected  to  agree  with  more  exact  calculations  or 
experiment,  although  any  approximations  so  obtained  would  probably  be 
consistent  with  linear-theory  estimates.  The  three-dimensional  profile- 
correction  factor  is  thus  seen  to  be  impractical  for  application  to 
arbitrary  profiles  because  (l)  limited  types  of  profile  -  plan-form  com¬ 
binations  were  treated  in  the  previous  three-dimensional  drag  analyses 
and  hence  there  are  insufficient  calculations  or  formulas,  or  both,  avail¬ 
able  and  (2)  drag  peaks  associated  with  linearized-theory  calculations 
would  lessen  the  accuracy  of  the  approximation  for  certain  ranges  of 
Mach  number. 


A  correction  factor  which  is  more  general  in  application  than  the 
three-dimensional. correction  factor  is  that  based  on  the  ratio  of  the 
drag  coefficients  for  two-dimensional  airfoils.  These  coefficients 
(based  on  Busemann  or  Ackeret  approximations)  are  readily  available  for 
many  airfoil  shapes  and  in  any  case  are  easily  obtainable  for  thin  profiles 
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of  arbitrary  shape.  (For  example,  see  formulas  in  section  V  of  refer¬ 
ence  13. )  Find  the  section  wave-drag  coefficient  c^  for  the  profile  s^ 

and  the  section  wave-drag  coefficient  for  the  profile 


following  profile-correction  factor 

(c*)s 


Sq;  then  form  the 


L(Ca>oJ 


Two-dimensional  flow 


An  expression  relating  the  drag  coefficient  of  the  desired  wing  C* 
to  known  coefficients  can  then  be  written  as  follows: 


NMl  "  N 


Ml 


(°d) 


N*o 


Two-dimensional  flow 


M 


(1) 


A- 


(t/c)x 


A-i 

Al 

(Vc), 


1 

(t/c)1 


where  the  subscripts  outside  the  parentheses  or  brackets  indicate,  as 
before,  the  wing  or  section  characteristics  and  Mach  number.  (it  is 
interesting  to  note  that  an  application  of  the  above  correction  factor 
for  some  wings  of  arrowhead  plan  form  with  biconvex  profile  may  be  found 
on  page  21  of  reference  l4. )  Of  course,  this  type  of  correction  factor 
cannot  be  expected  to  render  as  reliable  an  approximation  for  a  wing 
with  the  leading  edge  swept  behind  the  foremost  Mach  cone  (subsonic  leading 
edge)  as  for  a  wing  with  the  leading  edge  protruding  from  the  foremost 
Mach  cone  (supersonic  leading  edge). 

This  correction  f actor  was  used  to  estimate  the  drag  coefficient 
for  an  infinite  series  of  wings  with  the  following  characteristics: 

\  =  0.531*  A(3  =  4.65*  3  cot  A  =  1.375;  and  biconvex  profile  composed  of 
two  symmetrical  parabolic  arcs.  A  value  of  Cp^  for  a  family  of  wings 

that  was  geometrically  similar  but  had  rhombic  airfoil  sections  was 
obtained  from  reference  6.  Application  of  equation  (l)  yielded  the 
following  result:  Cp  =  4.88(t/c)2  cot  A.  This  value  compares  excellently 
with  the  theoretically  calculated  value  of  4.86(t/c)2  cot  A  (see  refer¬ 
ence  9).  This  remarkable  agreement  will,  of  course,  not  be  true  in 
general,  but  reasonable  estimates  should  be  obtained  for  most  cases  where 
there  are  no  fictitious  drag  peaks  associated  with  linearized- theory 
inaccuracies.  Fairly  good  approximations  may  especially  be  expected  at 
speeds  for  which  the  wing  leading  edge  is  supersonic. 
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In  this  connection,  some  discussion  of  the  results  obtained  in 
reference  10  would  be  appropriate.  Figure  14  of  that  paper  presents  the 
ratio  of  wave-drag  coefficients  for  two  wings  of  equal  thickness  ratio 
plotted  against  a  generalized  Mach  number  -  sweepback  parameter,  one 
wing  having  the  biconvex  profile  mentioned  in  the  previous  example  and 
the  other  a  rhombic  profile.  The  curves  of  figure  14  show  the  ratio, 

(based  on  three-dimensional  analyses  for  the  wings  considered)  to  vary 
from  approximately  1  to  1.5.  On  the  basis  of  the  profile-correction 
factor  proposed  in  the  present  paper,  the  estimated  ratio  would  be  1.33 
for  the  entire  Mach  number  -  sweepback  range.  The  estimated  value  for 
the  entire  range  thus  appears  to  be  fairly  satisfactory,  especially  when 
the  magnitude  of  the  wave-drag  coefficients  is  considered.  However,  a 
closer  examination  of  the  results  (see  fig.  12  of  reference  10)  will  show 
that  the  portions  of  the  curves  that  deviate  most  from  the  value  1.33  are 
largely  influenced  by  the  erroneous  drag  peaks  that  are  inherent  in  the 
linearized-theory  treatment  for  the  rhombic-profile  wing.  These  regions 
should,  of  course,  be  excluded  from  consideration  as  previously  discussed. 
Therefore,  the  estimated  factor  of  I.33  appears  even  more  satisfactory. 
Although  the  types  of  sections  considered  in  the  previous  examples  are 
restricted  to  biconvex  and  rhombic  profiles,  the  results  are  applicable 
to  a  considerable  variety  of  wing  plan  forms:  rectangular,  triangular, 
arrowhead,  sweptback,  and  sweptforward  wings  of  arbitrary  taper  ratio 
(conventional)  with  streamwise  tips.  Reasonable  estimates  of  the  drag 
are  obtained  even  for  wings  that  are  swept  well  behind  the  foremost 
Mach  cone,  although  no  generalization  can  be  made  regarding  this  agreement. 
However,  inasmuch  as  the  slope  of  the  airfoil  surface  must  be  small  every¬ 
where  it  appears  that  reasonable  estimates  of  the  wave-drag  coefficient 
should  also  be  expected  for  other  permissible  profiles  -  especially  at 
speeds  for  which  the  wing  leading  edge  is  supersonic. 

In  connection  with  the  previous  discussion  relating  to  the  erroneous 
drag  peaks,  it  is  quite  interesting  to  note  that  if  the  value  of  Cp-^ 

for  a  rhombic -profile  wing  is  obtained  by  using  the  value  of  Cp^  for 

the  parabolic-arc-profile  wing  at  a  speed  range  for  which  the  maximum¬ 
thickness  line  is  near-sonic  (Mach  lines  parallel  the  maximum-thickness 
ridge  line) ,  a  lower  estimate  for  Cp]_  is  obtained  than  is  calculated 
by  linearized  theory.  The  lower  value  would  agree  better  with  experi¬ 
mental  results.  This  agreement  is  due,  of  course,  to  the  fact  that  the 
effect  of  the  fictitious  drag  peak  associated  with  the  theoretical 
solution  for  the  rhombic-profile  wing  is  eliminated  by  using  the  correction 
factor  in  this  manner. 

In  addition  to  the  limitations  on  Mach  number  resulting  from  near- 
sonic-  or  sonic-edge  conditions  and  the  degree  of  applicability  for  wings 
with  subsonic  leading  edges,  the  usual  restrictions  of  linearized  three- 
dimensional  drag  analyses  are  applicable.  Regarding  the  geometric 
similarity  of  the  two  plan  forms,  the  sweepback  of  the  maximum-thickness 
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line  should  be  as  nearly  the  same  for  the  two  wings  as  available  calcu¬ 
lations  permit.  (The  leading-edge  sweepback  is  the  same  for  both  wings, 
as  indicated  in  equation  (l).)  The  proposed  correction  factor  should 
of  course,  be  applied  only  to  wings  with  section  profiles  that  satisfy 
the  thin-airfoil  restriction  that  the  slope  be  small  everywhere  on  the 
surface . 


CONCLUDING  REMARKS 


..  A  semiempirical  profile-correction  factor  is  discussed  which  enables 
the  estimation  of  the  wave  drag  due  to  thickness  at  supersonic  speeds 
for  three-dimensional  wings  with  arbitrary  airfoil  section  (subject  to 
the  restrictions  of  thin-airfoil  theory)  through  use  of  previously  calcu- 
a  ed  drag  coefficients.  It  is  felt  that  judicious  use  of  the  correction 
actor  should  yield  satisfactory  accuracy  for  many  plan-form  -  profile 
combinations,  especially  at  speeds  for  which  the  wing  leading  edge  is 
supersonic.  Application  of  the  proposed  correction  factor  to  experimentally 
determined  three-dimensional  drag  coefficients  presents  itself  as  an 
interesting  possibility. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va. ,  July  31,  1951 
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